Fluctuation-dissipation relations and energy landscape in an 
out-of-equilibrium strong glass-forming liquid 



Antonio Scala 1 , Chantal Valeriani 1 , Francesco Sciortino 1 ' 2 and Piero Tartaglia 1,2 
1 Dipartimento di Fisica e Istituto Nazionale per la Fisica della Materia, 
Universitd di Roma La Sapienza, P.le Aldo Moro 2, 1-00185, Roma, Italy 
2 INFM, Center for Statistical Mechanics and Complexity, Universitd di Roma La Sapienza 

(27 Sept 2002) 



Abstract 



We study the out-of-equilibrium dynamics following a temperature-jump in 
a model for a strong liquid, BKS-silica, and compare it with the well known 
case of fragile liquids. We calculate the fluctuation-dissipation relation, from 
which it is possible to estimate an effective temperature T e ff associated to the 
slow out-of-equilibrium structural degrees of freedom. We find the striking 
and unexplained result that, differently from the fragile liquid cases, T e ff is 
smaller than the bath temperature. 
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A central issue of statistical physics is to link the static properties of a thermodynamical 
system to its dynamics. A fundamental result in this direction has been the derivation of the 
fluctuation-dissipation theorem (FDT) [1], which links the response functions to autocorre- 
lation functions; however, derivations of the FDT require the system to be in equilibrium 
and the response to be in linear regime. 

Recently, a great deal of work has concentrated on the generalization of the FDT to 
out-of-equilibrium (OOE) situations [2]; in particular, it has been proposed to generalize 
FDT in the form [3] 



where A and B are conjugated observables, RA B {t,t') is the mutual response of the observ- 
able A to an impulsive conjugated field applied at t' , C AB (t, t') = (A(t)B(t')) - (A(t)) (B(t')) 
is the unperturbed correlation function of the observables and T is the temperature of the 
bath; X is supposed to be a function of Cab only and equals to one in equilibrium and 
linear regime. In recent analytical and numerical work, it has been observed that the out 
of equilibrium dynamics which follows a T-jump is characterized by an FDT relation with 
X = 1 at short times and X < 1 at long times [4]. It has been proposed to interpret 
the ratio T e ff = T/X as an effective temperature associated to the different timescales of 
the system [5]: while fast degrees of freedom are supposed to be thermalised so that T e ff 
equals the bath temperature T (and hence X = 1 and ordinary FDT holds), slow degrees of 
freedom are still reminiscent of the temperature of the system before being brought in OOE 
so that T e ff 7^ T (i.e. I ^ 1). In general, it is observed that for quenching experiments 
T e ff > T for the slow degrees of freedom; also in the case of a "reverse" quench, where 
the system is suddenly heated up, slow degrees of freedom "remember" the initial state 
being T e ff < T [6]. For fragile liquids, the results are consistent with the one step replica 
symmetry breaking (1RSB) scenario both for quenches [4,7] and for "crunches" [8] (sudden 
changes of volume); in this scenario, X becomes a step function equal to one at short times 
(T e ff = T) and equal to a constant less than one (T e ff > T) at long times. Exceptions to 
this behavior are pointed out in refs. [9-11]. 

The slow dynamics of supercooled systems have been considered also under a compli- 
mentary approach, based on the statistical properties of the free-energy [12] or the potential 
energy [13] landscape. In the landscape approach the free energy of a system is separated 
in configurational contributions, linked to the basins of attraction of local minima of the 
landscape, and in vibrational contributions linked to the number of states in such basins. 
It has been suggested that fast relaxations are related to the exploration of states inside 
a basin, while slow relaxations correspond to the exploration of different potential energy 
landscape basins. The topology of the landscape explored should therefore correlate to the 
dynamics of the system. Support of this interpretation for fragile liquids has been presented 
both for equilibrium [14] and out-of-equilibrium conditions [15]. 

Under the assumption that the basins explored in OOE are the same that would be ex- 
plored in equilibrium [16], it is possible to extend the landscape approach to OOE situations 
at the expenses of adding an additional parameter to the thermodynamic formalism. In 
OOE, the landscape-based liquid free energy can be written as [7] 



RAB(t, t') 



X(C A , B (t,t>))dC AB (t,t>) 
k B T dt> 



(1) 



(2) 
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where Ci S is the average energy of the local potential energy minima (the so called "inherent 
structures" [13]), the configurational entropy s con f counts the number of minima of energy 
Cis, fvib is the average free energy of a basin of depth e is , T is the temperature of the bath 
and T int is a Lagrange multiplier that extremizes /; in the equilibrium case, T int is just the 
temperature of the bath T. The extremum condition df/de is = gives an expression of the 
internal temperature T int 



that can be calculated from the equilibrium values of s con f and from analytical approxima- 
tions for f vib . In the case of a prototype fragile liquid, the binary mixture Lennard- Jones, 
it has been verified that in OOE conditions the configurational temperature T int calculated 
from the landscape approach and the effective temperature T e ff measured from the extension 
of the FDT coincide [7]. This coincidence is supported from the analytical predictions for 
the prototype fragile mean field system, the p-spin model, for which again T int = T e ff [17]. 

In the case of strong liquids, the situation is less developed. Only recently, a realistic 
model for liquid silica, the BKS model [18], has been studied in detail in equilibrium condi- 
tions. At variance with fragile systems, where the T dependence of the characteristic times is 
always super- Arrhenius, BKS-silica shows a high T region (T > 3330-fT) where the dynamics 
is super- Arrhenius and a low-T region (T < 3300/T) where the equilibrium dynamics recov- 
ers an Arrhenius behavior. For T > 3330.fr, the dynamics of BKS silica is well described 
by mode-coupling theory; at difference with other fragile liquids, the presence of local tetra- 
hedral structures and therefore of a strong oxygen-silicon-oxygen correlation, makes also 
three-point correlation functions to be taken in account in order to get quantitative agree- 
ment with the theory [19]. For T < 3330i^, BKS silica becomes strong, in agreement with 
the behavior of real silica. Interestingly enough, despite this strong-to-fragile crossover, the 
structural relaxation times of BKS-silica are well described in the whole temperature region 
by the Adam-Gibbs formula [20] as in the case of fragile liquids; moreover, free directions 
in configuration space vanish at the mode-coupling temperature [21]. Therefore, landscape 
theory is able to rationalize the equilibrium dynamics of fragile and strong systems despite 
the alleged diversity of their landscapes [22]. 

It is therefore important to check if T int = T e ff also in the case of strong liquids: if 
this would be the case, we could have afforded a simple thermodynamic understanding 
of OOE liquids in terms of two temperatures, the kinetic temperature T associated with 
the fast vibrational degrees of freedom and the internal temperature T int (or its dynamical 
counterpart T e ff) associated to the slowly evolving configurational degrees of freedom [23]. 

In our simulations we have employed the BKS-potential [18] with the parameters of 
ref. [20]. Simulations are for systems of A = 999 particles (N S i = 333 and N = 666) at 
fixed density p = 2.36g/cm 3 with a Nose '-Hoover thermostat. Long range forces have been 
calculated implementing the Ewald sums. We have quenched 280 independent configurations 
equilibrated at T = 3800 A" (well above the fragile-to-strong crossover at T = 3330 A" [24]) at 
two different final temperatures Ti = 2900 and T 2 = 2500; both these temperatures are in 
the region where BKS-silica has a strong behavior [24] . The Nose'-Hoover parameters have 
been selected to bring the kinetic energy of the system in equilibrium within 0.5 ps. We have 
measured two unperturbed correlation functions of the form {A(t + t w )B(t w )) , where A and 




(3) 
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B are conjugated observables with zero average in absence of a perturbation. Indicating with 
ff the coordinates of the /-th oxygen atom (/ = 1..No,Nq = 666), the choice A = p° coll = 
Ef° exp (ik ■ if) /sfN and B = (p° oll + 

Pcoii) corresponds to the dynamical structure factor 

Ccoii] the choice A = p° elf - = exp (ikr^ and B = (p° et f + p° & *eif) corresponds to the self part of 
the intermediate scattering function C S eif- In the case of the self correlation, averages over 
several ions have been performed. The total number of independent realizations (different 
configurations, different wave vectors, different ions) is ~ 30000 for the quenches at T x and 
~ 20000 for the quenches at T 2 , both for the self and for the collective case. 

Fig. 1 shows the evolution of the average inherent structure energy during the aging, 
following the T-jump and compares it with the known equilibrium values. The average ers 
decreases with time, while the system searches for deeper and deeper basins in the attempt 
to equilibrate. We note that the t{ S decrease does not depend on the bath temperatures 
selected, suggesting that the system is exploring configuration space with a saddle-dominated 
entropic dynamics [25]. For example, after 2 ps, independently from the bath temperature, 
the system is exploring basins which are typically explored in equilibrium at T ps 3650/T. 

By calculating the density of states of the explored basins and comparing it with the 
corresponding quantity evaluated in equilibrium, we confirm that the basins explored during 
the aging dynamics have the same density of states of the basins explored in equilibrium. 
Since in the case of BKS-silica, the density of states does not depends on ej S , in harmonic 
approximation df V ib/dei S = 0. This implies that, according to Eq. 3, the T int of the BKS- 
system when populating basins of average depth e is coincides with the temperature at which 
basins of depth e is are populated in equilibrium (Fig.l). It is therefore possible to calculate 
the internal temperature of the system inverting the equilibrium relations between ei S and 
T [7] (Fig. 1). At the waiting time t w = 2ps, we find that T int ps 36507^ for both the 
quenches, so that T int > Ti,T 2 and we expect to measure for the slow degrees of freedom a 
T e ff higher than the bath temperature. 

In order to measure T e ff, we evaluate the integrated response function x = 
V ft w dt'R(t,t'); from Eq. 1 the slope ^ of the parametric plot of x versus C yields 
Vo ^p = — Vo k B T ff - We have applied conjugate fields VoB0(t — t w ) of constant ampli- 
tude Vq after a waiting time t w = 2ps, significatively bigger than the characteristic time of 
the microscopic dynamics (see Fig. 2), and measured the integrated responses Xcoii = P° u 
and Xsdf = Pseif- We have chosen V = A.OkJ/mol for B = (p° coU + p° c l u ) and V = 5.0 J /mol 
for B = (p° se if + Pstif) in order to have a good signal-to-noise ratio in the response function 
while being still well inside the linear regime region [28]. We have averaged response and 
correlation functions over 300 wave vectors of length \k\ = 28.6nm _1 corresponding to the 
main peak of the static structure factor in order to have a good signal-to-noise ratio for the 
correlation functions [29]. We show in Fig. 2 the correlations and the response functions of 
the system. At the t w chosen, exploration of deeper basins has started as shown in Fig.l. 
This waiting time, which is much shorter than the infinite-time limit in which analytical 
predictions are derived, is consistent with the aging time chosen in most of previous studies 
of fragile liquids [30]. 

The parametric plots of x versus R for the quenches at the higher temperature T\ are 
unusual (Fig. 3, left panels), in the sense that although we know the system is aging, the 
slope of the plots is consistent with a T e ff = T along all the curves as in equilibrium FDT. 



4 



This behavior is similar to the results of ref. [10] for the Fredrickson-Andersen model, but 
contrasts with our expectation of a T e ff higher than the 7\. The parametric plots of x 
versus R at T 2 (Fig. 3, right panels) yields even more astonishing results: while at short 
times (high values of C) the slope ^ is equal to — Vojr^jr yielding the usual equilibrium 
FDT relation, in the long time region (low values of C), the parametric plot bends upward, 
i.e. T e ff is lower than the bath temperature. 

In conclusion, we have shown that for a realistic model of strong liquid, the aging dynam- 
ics is not consistent with the commonly observed 1RSB scenario. Moreover, the behavior 
of the fluctuation-dissipation ratio X contrasts with the possibility of interpreting T e ff as 
a temperature. From the point of view of landscape theory, this effect could be due to a 
profound and still to be understood difference in the way strong and fragile landscapes are 
explored. As for strong liquids barriers to relaxation grow slower than for fragile liquids 
while exploring deeper regions of the landscape, the possibility remains in fact open that ac- 
tivated processes are "too fast" and don't allow to describe the system in terms of long-lived 
quasi-states [31] and to define therefore a simple extension of landscape thermodynamics to 
OOE conditions. The results reported in this letter clearly contradict the general consensus 
built in recent years in the study of models of fragile liquids [4,7,8]. This pushes for a deeper 
theoretical understanding of the behavior of the strong glass-forming liquids. 

We acknowledge support from MURST PRIN 2000, INFM-PRA-HOP, INFM-Initiative 
Parallel Computing. 
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FIGURES 




FIG. 1. Left panel: equilibrium relation between ei s and T. Right panel: ej s as a function of 
t for the two quenches considered. The arrows show graphically the procedure that connects the 
value of ei s (t) to the equilibrium temperature r(e, s ) for t = t w = 2ps; landscape theory therefore 
predicts the temperature associated to the configurational degrees of freedom to be Ti nt ~ 3650 K 
in both quenches. 
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FIG. 2. Correlation and response functions for the collective (upper panels) and the self (lower 
panels) variables, at the two temperatures T\ = 2900 K (left panels) and T2 = 2500 K. Short time 
behavior corresponds to the fast vibrational dynamics inside a basin; the corresponding first decay 
of the correlation functions to a plateau is fast. Long time behavior corresponds to the exploration 
of new minima causing the slow decay of correlation functions from the plateaus. 
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FIG. 3. Parametric plots of the correlation and response functions for the collective (upper 
panels) and the self (lower panels) variables, at the two temperature T\ = 2900 K (left panels) 
and T2 = 2500 K; the slope ^ is proportional to —1/T e jf. Straight lines are the predictions of 
the equilibrium FDT theorem where the slope is proportional to — 1/T. At the lower temperature, 
T e f j < T at long times (low values of correlation) . 
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